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^—^ • Abstract 

> : 

en , We assume that the most important quark correlations are pairwise at 

I all baryon densities. We introduce correlated pairs by means of Bogoliubov 

' transformations which are functions of time and spatial gauge fields, in the 

, formalism of the transfer matrix with lattice regularization. The dependence 

T-H I on time and gauge fields allows us to enforce gauge invariance and other 

■ symmetries term by term in the transformed quantities. The resulting action 

• • . should be suitable for the description of multiquark mesons and baryons as 

. !_H I states of a quark and a diquark. We derive the quark contribution to the free 

^ ■ energy at finite chemical potential in a certain approximation. Its expression 

. cannot be evaluated analytically, but it has a definite sign. 



1 Introduction 



Our understanding of QCD, at least in the nonperturbative regime for the gauge 
couphng constant, strongly relies on numerical simulations which, indeed, have 
become more and more performing over the years. It would be desirable, however, 
to construct from first principles an approximate description suitable to capture 
efficiently what we believe are the essential features. And, possibly, introduce a 
perturbative scheme (in some parameter) in order to improve systematically the 
initial approximation. 

This kind of approach becomes more appealing when we wish to investigate 
QCD at finite chemical potential, because the numerical simulations of the fermion 
sector are plagued by the infamous sign problem. For example, in the fermion 
contribution to the free energy, this difficulty is revealed by the large accuracy 
needed in the evaluation of terms which (almost) cancel out among themselves q 

We suspect that the fermion contributions affected by sign fluctuations are 
due to states of high energy. In such a case, since these contributions to the 
free energy (almost) cancel out, an approximation which retains only the more 
stable fermion states, for any gauge field configuration, would reasonably solve the 
problem, simply because it neglects these fiuctuations altogether. 

There are several indications which might help in selecting such fermion states. 
At low baryon density, many authors think that diquarks are important substruc- 
tures in hadrons. Actually, in a historical perspective, baryonic constituents with 
diquark quantum numbers were already hypothesized by Gell-Mann [5], but as 
elementary constituents. Later, with various motivations, models of baryons con- 
structed in terms of one quark and one diquark have been investigated , and 
the diquark was regarded as a really composite state, even though, in practice, 
it was often treated as elementary. Subsequently, diquarks played an increas- 
ingly relevant role in the interpretation and explanation of several properties of 
mesons [H - IT8] . 

Also at very high baryon density diquarks appear of fundamental importance 
as the various phases of color superconductivity can be understood in terms of 
their condensation [T9H22] . 

A relatively smaller amount of work has been done at intermediate densities, 
where, however, it has been suggested that the structure of condensed diquarks 
might change with decreasing baryon density, while their size might shrink down 
to a dimension comparable with the average interquark distance. And this could 
explain the crossover from color superconductivity to Bose-Einstein condensation 
of molecular diquarks [23ti26] . 

^ For attempts to tackle this problem by analytic continuation in the chemical potential one 
can see [IHl]. 
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In conclusion, we deem that the results reported above H strongly suggest that 
two-quark correlations are very important at all baryon densities. In the following 
we shall take seriously this hint and we will assume that fermion states which are 
energetically stable, should always contain diquarks, in a condensed phase at high 
baryon density, accompanied by unpaired quarks at low baryon density. 

We must emphasize that when we talk about a diquark we mean only a pair 
of correlated quarks. Therefore, the simple presence of diquarks does not neces- 
sarily imply the existence of real or virtual bound states of two quarks, as the 
molecular diquarks of a Bose-Einstein condensate or the Cooper pairs of a color 
superconductive phase. Only above some critical values of the chemical potential 
such states can, eventually, appear. 

At a formal level, we will adopt the lattice regularisation and use the Kogut- 
Susskind formulation for fermions. 

As a first step, we need to identify positive and negative energy states of the 
Dirac lattice Hamiltonian. For this purpose we perform a first Bogoliubov trans- 
formation, which is equivalent to a Foldy-Wouthuysen transformation [27]. Bare 
fermions are replaced by quasiparticles in presence of a background field. This 
construction is really meaningful only for a class of gauge configurations that we 
shall call stationary meaning that the corresponding chromomagnetic field is in- 
dependent on time, and the chromoelectric field is vanishing [25]. We remark that 
these conditions are well suited to study also the effect of an intense background 
magnetic field on strong interactions, a problem considered of interest both at the 
level of the cosmological electroweak phase transition and for the heavy-ion col- 
lisions. And, indeed, numerical simulations have already been performed both in 
the quenched approximation (see, for example, [29]) and with dynamical fermions 
(see [5U] which contains also a detailed bibliography), in order to try to understand 
magnetic catalysis, i.e. the increase of chiral symmetry breaking induced by the 
magnetic background field. 

Afterwards, in this paper, we introduce diquarks by means of a second Bogoli- 
ubov transformation, so that diquarks will appear as Cooper pairs of quasiparticles. 

In the fermion contribution to the free energy, at fixed stationary gauge-field 
configuration, we can distinguish, in our formalism, a contribution Sbo from the 
vacuum, which we call bosonic, from the fermonic action of the remaining quasi- 
particles. At low temperature Sbo is dominating and becomes exactly the whole 
free energy at zero temperature. At vanishing chemical potential /i, it takes the 
form 

56o = -^trlng = -^5^1ng, (1.1) 



^Thc quoted references are only a sample of a vast literature, chosen to support our arguments. 
We apologize to Authors whose important contributions have not been acknowledged here. 
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where Lq is the size in the temporal direction of the lattice, and is therefore the 
inverse temperature, and Lo/2 is the number of blocks in which we have to divide 
the time-direction. The operator Q is related, as we shall see later, to the operator 
(see its explicit form in (13. 6p ) appearing in the definition of the transfer matrix 
of the quarks, and the Qi (respectively Ni) are its gauge-independent eigenval- 
ues corresponding to eigenstates that we enumerate by using the index i. More 
precisely 



For all states i these eigenvalues are, by inspection, real and larger than unity, that 
is InQj > 0. The result (11.11) coincides with what we got without the introduction 
of di quarks [28] . 

Let us now introduce a positive chemical potential. We find that Sf,o decreases 
according to: 



where the set P denotes quasiparticles states i such that InQj < 2/i. The above 
estimate has been derived by assuming that diquarks are formed by quasiparticles 
in P with a simple pairing structure, namely, for each i & P there is one and only 
one conjugate state p{i) G P, with Qp(i) = Qi. 

As the eigenvalues Qj's do not depend on the chemical potential, the fermonic 
number np, defined by 



exactly counts the number of paired quark states of the ensemble P. Thanks to 
the introduction of diquarks we have taken into account in the bosonic action the 
nonvanishing fermionic number. At increasing chemical potential the background 
field is depleted because of Pauli blocking. A short account of this analysis has 
already been presented in [31] . 

To derive (II. 3p we computed the bosonic contribution to the vacuum after the 
two Bogoliubov transformations which introduce, respectively, the background and 
the diquark field. If, on one hand, (II. 3p is recovered by the fields which satisfy 
a variational principle, the general expression, on the other hand, can be used to 
study multi-quark mesons and baryons as quark-diquark composites by taking into 
account the fluctuations of the background and diquark fields, along the lines of 
the expansion presented in p52j. 

For each fixed stationary gauge-field configuration we get from the condition 
log Qi < 2fi a sharp Fermi surface. After the integration on this space of gauge- 
field configurations we expect that the Fermi surface should be smoothed out. This 






(1.3) 
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is confirmed by a perturbative calculation in the gauge coupling-constant that we 
performed for large values of the chemical potential. In this case we get also a gap 
equation compatible with that obtained by standard methods [33] . 

We think it is useful to compare our results with the nonrelativistic ones. From 
the technical point of view our formalism is a fermion number conserving extension 
of the theory of superconductivity developed by Bogoliubov and Valatin [3^IB5] 
which violates this symmetry. The enforcement of fermion conservation in many- 
body theories can indeed be achieved by allowing time-dependence of the Bogoli- 
bov transformation [36]. In the saddle point approximation, however, one gets 
a formulation close to the quasi-chemical equilibrium theory of superconductors 
developed by the Sydney group [HZ], in which fermion number is explicitly pre- 
served. Since the latter approach is in our opinion more transparent than that 
of BCS and Bogoliubov- Valatin from a physical point of view, establishing a con- 
nection between superconductivity and superfluidity, we report a brief account of 
both methods in Appendix A. 

The plan of the work is as follows. In Section 2 we report some considerations 
about the fermion determinant with special regard to its form in the presence of 
stationary gauge fields. In Section 3 we report the definitions and notations we 
will use and in Section 4 the time-dependent gauge field-dependent Bogoliubov 
transformations. In Section 5 we write and solve the saddle point equations, in 
Section 6 we perform the perturbative expansion in the gauge coupling constant 
for large chemical potential and we conclude with some remarks in Section 7. 

2 The sign problem 

In this section we will review the well known sign problem which affects numer- 
ical simulations in presence of fermions in many problems and in the particular 
case which is interesting for us, namely QCD at finite density. 

We begin from the expression of the grand-canonical partition function of QCD 
directly in continuous space-time. Formally it can be represented as a path integral 
in euclidean space 



where A represents the gauge fields, Sg is their pure action, and the fermion 
partition function is given by a Berezin integral. 



where the fermion action S'i;'[A] is bilinear in the Grassmann variables 







5 



The explicit integration on the fermion fields provides the so-called fermion 
determinant 

Zf[A] = det[f + m + /i7o] (2.3) 

where ^ is the contraction of the covariant derivative, which depends on the gauge 
fields, with the Dirac 7-matrices, m is the fermion mass and fi is the chemical 
potential. 

Remark that as 

{f + m + fx^oY = - f + m + fj.*jo (2.4) 
75(y + m + /i7o)75 = - y + m - ;U7o (2.5) 

the fermion partition function is necessarily real both when /i is vanishing or purely 
imaginary, because under these conditions 

Zf[A] = det 75[y + m + fi^oH = det[f + m + /i7o]t = Z*p[A] • (2.6) 

In the Weyl (chiral) representation for 7-matrices 

where a are the Pauli matrices, the relevant matrix takes a partitioned form, 

y + m + /i7o= ,--v7 ' (2-8) 

which is particularly suitable to reduce the evaluation of the determinant to a 
space of half dimension, indeed, 

Zf[A] = det m'^ - {Vo + fi + ia ■V)(yo + fx-ia -V) . (2.9) 

It soon follows that, in the case of fi vanishing or purely imaginary, not only the 
determinant is real, but it is also non-negative. Indeed, if we set 

X = iVo + ifi + a -V , (2.10) 

then 

Zf[A] = det(m2 + X^X) > detX^X > 0, (2.11) 

where the first inequality becomes an equality for vanishing mass and the second 
whenever X has a vanishing eigenvalue. 

But for real, non-vanishing, chemical potential the fermion determinant ap- 
pears, in general, to be complex. A more detailed analysis can show that it is 
possible to combine gauge configurations in pairs so that the sum of the determi- 
nants be real [38]. It remains, however, a possible sign, which is the problem in 
dynamic Monte Carlo simulations where fermions are integrated out, because the 
resulting factor cannot be used as a positive weight which drives the importance 
sampling of gauge configurations. 
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2.1 More on the fermion determinant 

In the Introduction, we have already seen that we are interested in station- 
ary gauge-field configurations in which the chromoelectric field vanishes. In this 
Subsection we will investigate what are the consequences for the evaluation of the 
fermion determinant when we restrict ourselves to such gauge configurations, that 
is such that 

[Vo,V] = (2.12) 

the spatial covariant derivatives commute with the temporal one. Within this 
ensemble of configurations, the fermion determinant becomes 

Zf[A] =det[m2-(a-V)2-(Vo + /i)'] = det[H^ - {Vq + fif] (2.13) 

where the square of the Hamiltonian H, which is 

H^^m''-{a- Vf , (2.14) 

depends only on the spatial components A of the gauge field and is positive definite. 

The operator (Vo + /i), instead, depends only on the temporal component Aq 
of the gauge field. It is quite simple to see that if it has, for real a complex 
eigenvalue A, it has also the conjugate one, and therefore the fermion determinant 
can be written as a product of positive terms 

MA] =11 n {h'-X'){h'-X')>0, (2.15) 

h A:S(A)>0 

where /i^ are the eigenvalues of H^. 
A few remarks: 

• the determinant becomes a function of /i^, that is it does not depend on the 
sign of /i; 

• the determinant does not feel the sign of Aq, so that the average values of 
the Polyakov lines in positive and negative time directions are the same; 

• in the absence of the chromoelectric field, the temporal component Aq of 
the gauge field, which is responsible for the implementation of the Gauss 
constraint, appears only in the fermion determinant. The previous discussion 
implies that 

J [dAo] Zf[A] > . (2.16) 
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When using lattice regularization, as we shall do in the following of this paper, 
strictly speaking, condition f l2.12p cannot be used. Covariant derivatives must be 
replaced by their discretized versions, as for example, defined according to 

V+0(x) := f/^(x)0(x + fi)- 0(x) (2.17) 
V;</)(x) := <j){x) - Ulix - fi)<j){x - fi) . (2.18) 

The curvature fields are computed by taking the commutators of the covariant 
derivatives, so that, for example 

[V+, V+] <P{x) = [U,{x)U,{x + fi)Ul{x + u)Ul{x) - 1] (2.19) 

U^{x)U^{x + v)<P{x + fi + u) (2.20) 

= [U^^,{x) - 1] U,{x)U^{x + z>)0(x + fi + v) (2.21) 

[V^ , V+] <P{x) = [1 - Ul{x - fi)U,{x - fi)U,{x - /i + u)Ul{x)] (2.22) 

U^{x)Ul{x - fi + v)(l){x - fi + v) (2.23) 

= [1 - f/_^,.(x)] U,{x)UUx - /i + i))0(x - A + i>) (2.24) 



where U±^^±^{x) are the four plaquettes open at the vertex x in the plane (/x, z/). 
They transform, under a gauge transformation, according to 

U±^,±u{x) ^ n\x) t/±^,±.(x)^](x) . (2.25) 

The vanishing of the chromoelectric field means that the spatial link configurations 
at successive times are related by 

Ui{x + 0) = Ul{x)U^{x)Uo{x + i) (2.26) 

so that the open spatial plaquettes, which determine the chromomagnetic fields, 
are constrained to be gauge equivalent, indeed 

U±i,±jix + 6) = f/c|(x) f/±i,±,(x)?7o(x) , (2.27) 

and therefore their trace is time invariant and they give one and the same contri- 
bution to the gauge-field action at all times. 

3 Definitions and notations 

The partition function of lattice QCD can be written as 

Z = j [dU] expi-SclU]) Zf[U] , (3.1) 



where [dU] is the Haar measure over the Wilson hnk variables U, that lives in the 
gauge group, Sq is the Wilson action for the gauge fields and Zf[U] is the fermion 
determinant. For our purposes we shall make use of the operator formulation, so 
that 

Lo/2-l 
t=0 

In the above equation Lq is the size extension of the lattice in the temporal direc- 
tion, so that it is also the inverse temperature, T is the fermion transfer matrix 
which acts in the Fock space of fermions and Tr^ is its trace . We shall make use 
of the Kogut-Susskind formulation for lattice fermions, so that fermion fields live 
on blocks of size twice the lattice spacing. The index t labels the blocks along the 
temporal direction. 

The expression of the transfer matrix in the gauge f/o = 1, but for Wilson 
fermions and in the particular case r = 1 for the Wilson parameter, was given by 
Luscher [39] , who proved also its positivity. See also [IQIIIT] for the generalization 
also to different values of the parameter r. The extention to Kogut-Susskind 
fermions, in the so-called spin basis, was given in [321113] • We shall use, instead, the 
flavour basis because a simpler transfer matrix is avalaible for this formulation [H] . 

Without fixing the gauge, the transfer matrix, at nonzero chemical potential 
/i, can be written as 

Tt,t+i ■= f} Vt exp(2 /i h) ft+i . (3.3) 
where n is the fermion number operator 

h := M^-u — -O^-O , (3.4) 

(the sum on all the indices is understood) with and v\ creation operators of 
fermions and antifermions, obeying canonical anti-commutation relations and 

Tt = ex-p[vNtu\ , Vt = exp[uHnU(i^tU + InU^^v] . (3.5) 

The matrices Nt are functions of the spatial link variables at time t. More precisely 

N = -2 (70 ® 1) |m + ^(7, ® 1) + Pj^VJ"^] | (3.6) 



where 

) = ^ ( 1 ® 1 ± 7,75 ® htj ) (3.7) 
are projection operators, 7^ and are Dirac and taste matrices. 
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are covariant derivatives, Tj arc forward / backward translation operators of one 

block of size twice the lattice spacing and Uj the j-th component of U, the spatial 
link variables associated to the blocks. 
The operators 

P± = ^(10 1^7075 ®t5to) (3.9) 

project on the components of the fermion field which propagate forward or back- 
ward in time 

u = P+V 

= P_ip. (3.10) 

The symbol "tr" denotes the trace over fermion-antifermion internal quantum 
numbers and spatial coordinates (but not over time). We introduce the notation, 
which we will use for any matrix A 

tr±A := tr (P±A) . (3.11) 

Finally we will denote by Tq^^ the forward and backward translation operators of 
one block, that is two lattice spacing, in the time direction 

[To'^\i,t2 = ^t2,h±i ■ (3-12) 



4 Time-dependent Bogoliubov transformations 

We evaluate the trace of the fermion transfer matrix in a basis obtained by 
performing Bogoliubov transformations on the coherent states 

\a,p) ^exp{-av) - pv^)\0), (4.1) 

where the a,/3 are Grassmann fields. 

In a first transformation, we introduce quasiparticlcs operators a, (3 which have 
the same fermion number as the original operators il, v 

a = R^it - J^^ v^) , /5 = (v + ^t jrt) (4 2) 

where 

p = (l + J^t^)-i ^=(i + ^^t)-i. (43) 

The uppcrscript circle denotes the involution defined by the above equations. The 
new operators satisfy canonical commutation relations for any choice of the matrix 
J^. The vacuum of the new operators is a condensate of this composite boson 

I J-) = exp(^t) |o) (4.4) 
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where 

By a second Bogoliubov transformation, we introduce new quasiparticle oper- 
ators a which have the same fermion number as u 

a = r^ [a- V'^a^) , (4.6) 

where 

r = (l + pt2))-i (4.7) 

and the bosonic field represented by the antisymmetric matrix T> has fermion 
number two. The corresponding operator 

pi" = dTpt^t (43) 

will represent diquarks. 

The vacuum of the new operators is 

|P, J^) = exp Q^"^) exp(J^T) |o) = exp (^a^V'^a^ exp [u^T^'^) |0) , (4.9) 

namely a condensate of Cooper pairs of quasiparticles in a background field. 
If we perform a gauge transformation 

^(x) ^'{x) = g{x) tp{x) (4.10) 

both components u and t)^ transform in the same way, that is 

u{x) —J- u'{x) = g{x) u{x) , v\x) — v'\x) = g{x) v\x) . (4-11) 

In order to get that also a transforms in the same way we need, because of (14. 2p . 
that the matrix appearing in the first transformation at time t transforms according 
to 

( J-;).,y ^ ( J-f ),,y = g{t, x) ( J-/).,y g\t, y) . (4.12) 

As a consequence, JF and the states \ J^) are gauge invariant. 

If we now demand that also a transforms as a we need, because of (14. 6p . that 
the matrix appearing in the second transformation at time t transforms according 
to 

(Dj).,y ^ {V'X,y = git, x) {vXy g{t, y) , (4.13) 

which also implies that V is gauge invariant. Therefore under the conditions (14.12^ 
and (I4.13P all the states J^) are gauge invariant. 
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The trace of the fermion transfer matrix, after the first transformation, was 
represented in [271[28] as a Berezin integral, with the result 

Zf = /[rfac/a*rf/3rf/3*]e-^™(^)-^'*'("'^;^) (4.14) 

where the Grassmann variables satisfy antiperiodic boundary condi- 

tions in time. In the above equation S'qp(a, /3; J-"), which is the action of quasi- 
particles, takes the form 

Lo/2-l 



+ a\{^t - Ut)at^x - A+i(Vt - Ut) /3;J (4.15) 
where the covariant derivatives are defined as 

V. := \ [f^ Uu - , Vt := ^ (e-^^ < - T^) , (4.16) 

and Tq^'* are translation operators of one time block defined in (13.121) . The presence 
of the factors 2 is related to the fact that neighboring blocks stay at two lattice 
spacings. 

The explicit expressions for the mesonic action Sme{^)i for the Hamiltonians, 

o 

respectively of the fermions and the antifermions, T-Lt and l-Lti and for the mix- 
ing terms between quasi-particles and quasi-antiparticles /j+i'' and Z^^'^'' will be 
reported in the next Section. 

In the present work we are interested in the trace of the fermion transfer matrix 
after the second transformation. At this stage, as we will concentrate on the leading 
contribution in the nilpotency expansion, we restrict ourselves to the bosonic part 
of the action 

Sbo = Sme + Sdq , (4-17) 

and we disregard the action of the new quasiparticles appearing after the second 
transformation. More precisely we shall restrict to 

/ c ^ TT {^t,^t\Tt,t+i\T^t+i,^t+i) //I io\ 

exp(-S.):= n (V^HVuT,) ' 



t=0 



Srae) which cau be obtained by putting = at all times t, is what we got 
previously, and the extra term Sdq-, that we call the diquark action, is derived in 
Appendix C, and takes the form 

^ Lo/2-l 

Sdq = i^ J2 + ^*^0 - In (l + '^''^t Qt+i,t AVi Qt+i,) ] (4-19) 



i=0 
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where the matrix Qi^i ^ is defined by 

QT+iy-=Uo,t-2e-''^nt, (4.20) 

and we denote by ^ the transpose of the inverse of the matrix Qt+i,t (which 
is also the inverse of the transpose). 

In the following we will work in the gauge f/o,i = 1, for all times t, and we 
should impose the Gauss constraint on the states \Vt,J^t)- But we don't need 
it, because, as a consequence of their gauge invariance, they are left invariant by 
the Gauss projector Vg, which takes into account also the transformations fl4.12p 
and dUSD, that is 

VG\'Dt,J't) = \'DuJ't) (4.21) 

for all times t. 

It is important to call the attention of the reader on a relevant difference 
between the conditions fl4.12p and f l4.13p for the gauge invariance of, respectively, 
J-' and T>. Indeed, in order to ensure the gauge invariance of the operator j^, which 
acts on the fermion Fock space at time t, it is enough to allow a dependence of J-" 
from the spatial-link variables f/k,t, such that, under a gauge transformation 

with 

iJ'I[UL,t])^,y = iJ'l[9tU^,t9l]).,y = 9it,^) (J-/[f/M])x,y^?^(t,y) . (4.23) 

This indicates that in the matrix element the dependence on the link variables is 
by strings Ff (Cx,y) of products of link variables at time t along paths Cx,y between 
the positions x and y, which realize the parallel transport and transform according 
to 

TtiC^,y) -> git, x) Fi(Cx,y) g\t, y) . (4.24) 

An attempt to proceed in the same way with the operator T> does not work. 
For example, for the case in which the gauge group is SU (3), in the simple proposal 

[Vlr^'^ = e,,y, (Ff '^(C^,x)rf ^(C^,y) - Ff ''(Cw,y)Ff «(C^,x)) d%t, w) , (4.25) 

where D is explicitly antisymmetric, as it must be, and we have also indicated 
the gauge-group indices, the gauge-invariance condition (14.130 forces d'^it, w) to 
transform, under gauge transformations, as a quark field. Therefore, in order to 
implement the gauge-invariance of P we must introduce a new dynamical field for 
diquarks, which acts as a compensating field, according to the general discussion 
on symmetry breaking we presented in PH]- 
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4.1 Expressions for the mesonic and quasi-particles action 



By using the matrix J-", which defines the first Bogohubov transformation, and 
the matrix N, which appears in the definition of the transfer matrix (13. Sp . we 
introduce the shorthands 



J-jv, t := 1 + N^J^t , , := 1 + J^^n} 

in terms of which we define the expressions 



Et+l,t ■—J^N,t+lUo,tJ^N,t + ^t+l^(),t^ ^ 
o o /o \to /o\t 

Et+l,t '■= J^N,t Uo^t ( J^N,t+l I + Tt Uo^t ( 



(4.26) 

(4.27) 
(4.28) 



With the help of the definitions of R and R, given in (14. 3p . we report now the 
expression for the Hamiltonians for the fermions and the antifermions 



:= I e'^ (Uo,t - R^ E-_,\,R;'\ 

m := - e-'^ (Ul- R,_,\ E;\, R, ' 
and for the mesonic action 

io/2-l Lo/2-1 

S^e{J^):=- J2 tr+ln(i?jEt+i,t) = - ^ tr+lnQ 



(4.29) 
(4.30) 

(4.31) 



i=0 



t=0 



where we used the definition (I4.20p . and for the the terms which mix quasi-particles 
with quasi-antiparticles 



r(2,l) ._ 1 ° 5 ^° 

h - 2 Rt 



Rt— EtJ_i J^NA-l 



Rt — Uo^t [J^N,t+i) Et+i 



o 1 

Rl 



(4.32) 
(4.33) 



As a consequence of the definitions (I4.20p and (14.291) we get the relation 

Qt,t+i ■■= Ri,Et+i,tR! ■ (4.34) 



5 Saddle point equations 

We determine the values of the background and diquark fields by minimizing 
the bosonic part of the effective action. The stationarity equations are 

d d 

Sbo = T;zr- Sbo = , (5.1) 
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We remark that in the nilpotency expansion we have to integrate, on an arbi- 
trary measure, both on the background and the diquark fields. Then stationarity 
equations become saddle point equations in the nilpotency expansion. There is a 
subtlety. In the latter case we should also add the equations 

d d 

" Sbo = 0, (5.2) 

and consider the fields and their Hermitian conjugates independent in the varia- 
tions. It will turn out that the solutions are Hermitian fields, so it will be enough 
to consider the solutions of (15. ip by keeping fixed the Hermitian conjugate fields. 



5.1 In the absence of the diquark field 

In the absence of the diquark field, that is whenever = for every time t, 
the background field has already been determined [28]. We outline its derivation. 
We look for solutions stationary in time, as appropriate to the vacuum. If J-" is 
stationary, the elementary bosonic fields coupled to the fermions which enter its 
expression should also be stationary. In gauge theories J-" must certainly depend 
on spatial link variables ?7k(^,x). Stationarity in time for gauge fields can be 
formulated in a gauge covariant way by requiring that these fields evolve according 
to gauge transformations, that is 

f/k(t,x) = l^,|,f/k(0,x)iy,^,^t. (5.3) 

As a consequence, the chromomagnetic contribution to the pure gauge-field action, 
namely, the trace of spatial plaquettes, does not depend on time. 

Accordingly, the matrices J^t and Nt are related to that at time t = 0, that is, 
if Jo = J- and Nq = N then 

jr^ = w}7Wt , Nt = W^NWt . (5.4) 

We still wish to set the contribution of the chromoelectric field to the gauge-field 
action, namely, the trace of spatio-temporal plaquettes, to be independent on time. 
We have been able to arrive at a stationary solution for J-" only with the particular 
choice 

Wt+i,^ = f/o(0, x)f/o(l, x) . . . f/o(t, x) (5.5) 

which lets the contribution from the chromoelectric field vanish at all times. 

Under these conditions the saddle point equations for the background field 
become independent of time 

jr = AT + jr (jr^)-i . (5.6) 
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The relevant extremal solution for the background field is 



= A^(2A^^A^)-^ N + ^/{Nmyl^Wm . (5.7) 
This is also the solution of the equations [27] 

This means that at the minimum of the vacuum energy there is no quasiparticle- 
antiquasiparticle mixing. (This is in close analogy to the case of the Bogoliubov 
transformation in the BCS theory, as explained in Appendix A. 2. But needless 
to say, unlike the latter these terms do not violate any symmetry). Hence at the 
saddle point, the effect of the Bogoliubov transformations (14. 2 p is analogous to that 
of the Foldy-Wouthuysen transformations which separate positive from negative 
energy states in the Dirac Hamiltonian [27 \ . 

The time evolution of the quasiparticle Hamiltonians is slightly different: 

nt = w}HWt+i , ■Ht= wl, n Wt (5.9) 

and similarly 

QT'i,t = W^Q-'Wt+i- (5.10) 
At the saddle point, the quasiparticle and antiquasiparticle Hamiltonians at time 

o 

t = 0, respectively H and "H, are simply related by 

o 

2e-2^W = 2e2^?7= 1 -g~^ = 1 - (5.11) 
They are Hermitian functions of N'^N, and the vacuum energy is 

5^e = -Ytr+lnQ. (5.12) 

5.2 In the presence of the diquark field 

Now we rewrite the diquark field action exploiting the time dependence (15. 3p 
of the spatial-link variables, but we will not need the explicit expression of the 
background field. So we will write, in addition to (15.101) where Q = Qi^, 

Vt = WfVWt (5.13) 

where V = Vq. Then the time dependence disappears from the bosonic action 



5bo = ^tr+|-lng + iln(l+Dtp) _ 1 In (l + e^^ P Q-^ 1?+ Q-^) | . 



(5.14) 
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We first remark that all the dependence on the background field T is now contained 
in the dependence from Q, so that 



d7 dQ d7 

As the equation 



0. (5.15) 



does not involve the diquark field T> and its relevant solution for T is exactly the 
T given in (15 .Tp . the background field does not depend on the diquark one. 
The stationarity equation for the diquark field is 

V = e''^Q-^VQ-^ . (5.17) 

As the bosonic action is gauge invariant, but the field T> is not, the stationarity 
equation determines only the class of T> equivalent under gauge transformations. 

In order to analyze the solutions of (I5.17P we construct the diquark structure 
function in the basis of eigenstates of the quasiparticle Hamiltonian for given gauge 
field configurations, which according to (15.111) are also eigenstates of Q 

Q\i) = Qi\i). (5.18) 

Remark that, although the operator Q is nonlocal, its eigenvalues are simply re- 
lated to those of the local operator A^, as shown in (11.21) . This explicit relation 
shows also that they are all real and greater than unity. 
The saddle point equations then become 

Vij = e^^Q-^VijQ-^\ (5.19) 

First we notice that this equation can only determine because any possible 
phase factor cancels from both sides and we can restrict to the case in which 
Pjj is a nonnegative real number. Apart from the trivial solution Dij = 0, the 
saddle point equations are satisfied for arbitrary Vij provided Qi = e'^^Qj^, in 
which case the contribution to the action vanishes. The relevant minima of the 
action are then reached on the boundary of the range of the \T>ij\, namely when 
these are zero or infinity. In the first case the Bogoliubov transformationtion 
is the identical transformation, in the second case it interchanges creation with 
annihilation operators. 

So we must determine in which way to chose between zero and infinity for any 
given eigenstate of Q. We will not search for the most general form of the P. 
Instead we restrict ourselves to the following canonical form 

'^ij = ^j,p{i) ^« ' = '^p{i) (5.20) 
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in which any quasiparticle state i is associated to one and only one conjugate 
state p{i). In the many-body language this is called simple pairing. Notice that, 
since the Vij are antisymmetric, they can always be cast into this form, but at the 
price of a unitary transformation which would change the form of the quasiparticle 
hamiltonian. For the same reason, we observe that, since the are functions of 
the spatial link variables, pairing can be simple only in a given gauge. 

We will denote by P the set of the states i for which l^jl = oo. It follows that 
if i e P also p{i) e P. Therefore 



^ - 5; In + 1 5; [in \V,\^ - In (e^- m^Q-^Q-^, 

= ^ |-^lng,-5^[2/x-lnQ,]l (5.21) 

assuming = Qp{i), at least when i & P. For given chemical potential this 
action is minimal if Qi < e^^ for each state i & P. The contribution originating 
from diquarks in the states i E P cancels the contribution to the action Sf,o of 
the corresponding components of the background field. This has a simple physical 
interpretation: because of Pauli blocking the states occupied by diquarks are not 
accessible to fermions in the background field. At/j, — OasQ>l the set P is 
empty, so that Vi = for any i. Increasing the number of fermionic states in 
the background field is progressively depleted, until a deconfining transition will 
take place. 

We can perhaps understand better this result if we rewrite the bosonic energy 
in the form 

I i i&P ) 

Taking the formal continuum limit we get that Hi < /j, for the states i e P and 

Zf ^ e^(^'^'"^*e^(^'-'')) . (5.23) 

We remind that the eigenvalues Hi are functions of the gauge fields. Integrat- 
ing over the gauge fields with the pure gauge fields weight will smooth out the 
distribution of the values of iVi] — 0, 00. 



5.3 Variational character of the saddle point approxima- 
tion 

It is important to notice that the expression of the bosonic action can be 
obtained by evaluating the partition function in a Fock space which contains only 
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the state J-"). At zero baryon density, P = 0, a variational calculation was 
performed in [32] and, subsequently, it was recognised to provide the vacuum 
contribution obtained by a suitable Bogoliubov transformation [15] . Here we sketch 
the proof of this property in the more general case. 

Let us start from the expression (13.21) for the fermion determinant and choose 
the gauge Uo^t = 1 for all times t, imposing the Gauss constraint on the states by 
means of the Gauss projector Vg- Under stationarity conditions for the gauge-field 
configurations we get 



/ 



r 



Uo,, 



(5.24) 



where we don't need to remember the time indices for the transfer matrix. 

Our transfer matrix T at Uo,t = 1 is positive definite, indeed T ~ T^T, and we 
simply get the inequality (I2.16P on the lattice 



[dUo] Zf>0 



(5.25) 



Thanks to the positivity of the transfer matrix we shall now use the inequality 



Tr^ Vg T 



> 



Uo,t=l 



{A\A) 



(5.26) 



valid for any gauge invariant state |^) in the fermion Fock space. By choosing 
1^) = [D,T) we get 



[dUo] Zf > 



2 



= exp [-Sboi'D,!')] > exp [-SmdJ^)] 

(5.27) 

which shows that the fermion determinant is approximated from below by the 
exponential of minus the bosonic action S^o, which is always greater than the 
exponential of minus the mesonic action, obtained by putting V = in Sbo- 



6 Perturbative expansion in the gauge coupling 
constant 

At sufficiently high baryon density an expansion in the gauge coupling constant 
can be justified. We proceed as if a nonvanishing background field will survive. 
The study of the densities at which the background field vanishes is deferred to 
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a future work. This Section can then be considered propedeutic to such future 
work. But it will also give us the possibility of a critical discussion of perturbative 
calculations in the gauge coupling constant and of a comparison with the non 
relativistic case. 

Under our assumptions the diquark effective action does not depend on the 
temporal links, but it depends on the spatial ones through the matrix Q and the 
diquark structure functions V. In order to determine T> variationally, we expand 
Q in powers of the the gauge coupling constant 

Q-^^l + A + gB + g^C. (6.1) 

Setting 

V^V 1 , , 

and expanding the In with respect to lattice spacing and gauge coupling constant 
we get 



S,o^ - '^tii^nQ + p{A + g B + g^ C) - ]^g''pBpB + ]^ g^^B^j^ B | 



(6.3) 



Notice that ip is antisymmetric. The above expression is essentially identical to 
that of the many-body theory reported in Appendix A, Eq. flA.7p . in which the 
first term is the vacuum energy (independent on the chemical potential), the sec- 
ond the kinetic energy, the third the density-density interaction and the last one 
the interaction between fermions in one and the same Cooper pair. Such close 
correspondence exists because we constructed the diquark field in terms of quasi- 
particles. 

At this point we assume simple pairing, so that 

{V^V).. = 6ij\Vi\^ , Pij = 6ijPi , {ip)ij = 5j^p{i) eij (6.4) 
where ^ 

We remind that we use the label i for all quasiparticle labels, that is the position 
vector X, the Dirac a, the favour / and the colour a indices. As a concrete example 
of simple pairing for the case of 2 flavors we can assume 

1^y.,aji,a,y,l3f2,b = e/1/2 ^3ab ^a/B ^x.y^x (6.6) 

x,a,/i,y,/3,/2 (4)y (6.7) 
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similarly for C (6-8) 

where k labels the spatial components and the A/ are the gauge group generators. 
The effective action now becomes 

5*60 ~ - Y '^\^'^Qi + Pi{^ + 9 B + C)ii - ^ g'^pi BijpjBji 

i ^ 

+ 1(^*A, + A*^,)} (6.9) 

where ^ 

= - g'^(^i,p{i) E ^k,p{k) Bik Bp(i)p(k) i'k (6.10) 

k 

is the celebrated gap function. Usually in many-body problems the density- density 
interaction is small but complicates the variational equation, and for this reason it 
is accounted for renormalizing phenomenologically the chemical potential We 
will adopt this criterion. Then we will integrate over the gauge fields with the pure 
gauge measure, and we will denote by (...) such average. Note that in the expansion 
we have disregarded the dependence of the diquark structure functions on the 
gauge fields. This can be justified only if the diquark is approximately pointlike, 
namely if its mean square radius is much smaller that the average interquark 
distance (see ref. [2D]). We will say a little more about this in the last Section. 
Bearing in mind that {B) = because B is linear in the gauge fields, we define 

Ki = Pes + iA + g^B%i -IYI 9^Pk{B^k B^i) . (6.11) 



2 

k 



The averaged effective action is 



{s,o) ^ - Y E + - /^^ff) + \ + ^r(A.))| . (6.12) 

Variation with respect to V (omitting the symbol of average on A)gives 

2{K, - pes)V, + ^ {Vj - 1) A, = (6.13) 

which has the solutions 

^i = ± — (6.14) 
2v/(i^,-/i)2 + |A,|2 

Inserting them into the definition of the gap function fl6.10p we get the gap equation 

2 "-^ {Kk - pY + \^kr 
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This expression has two specific features which agree with the standard result [T9] - 
[22]. The first is that in both cases the dominant contribution to the supercon- 
ducting gap comes from chromomagnetic fields. The second is that these fields are 
quasistatic in the standard approach, while in the present one they are completely 
time independent. 

7 Summary and future perspective 

We have investigated QCD guided by the theoretical indications that two 
quarks correlations are important at all baryon densities. We introduced such 
correlated pairs by means of Bogoliubov transformations in the formalism of the 
transfer matrix with lattice regularization. We performed a first transformation 
which produces a background field and quasiparticles with the quark quantum 
numbers. A second transformation yields the diquark field in terms of quasiparti- 
cles. Evaluating the trace in the partition function using a fermionic basis obtained 
by these Bogoliubov transformations on fermionic coherent states we got an effec- 
tive action of the system exactly equivalent to the original one. At variance with 
previous use of the Bogoliubov transformations we let them to depend on time 
and on spatial-link variables. This makes it possible to enforce for quasiparticles 
the same symmetry transformations as for quarks, and thus perform the saddle 
point approximation keeping gauge invariance manifest. The construction of the 
diquark field in terms of quasiparticles constitutes another innovation. 

We have investigated the effective theory at the zeroth order of a nilpotency 
expansion, namely an expansion in the inverse of the number of fermionic states 
in the structure functions of the composites, called the index of nilpotency. This 
means that we derived the effective action in a saddle point approximation, that we 
have shown to be equivalent to a variational calculation, in which the free energy 
is minimised with respect to background and diquark fields. 

The solution for the background field in the absence of diquarks was found in 
previous papers. According to this solution the QCD vacuum is a dual supercon- 
ductor in which the chromoelectric field is totally expelled from the vacuum and 
the fermion Fock space contains quasiparticles only in the form of point-like color 
singlets [28] . 

In the present work we have solved the saddle point equations in the presence 
of a diquark field. In order to describe multiquark mesons and baryons as bound 
states of a quark and a diquark we should include the fiuctuations of this field. We 
have restricted ourselves, however, to the study of diquark condensation at finite 
chemical potential. At fixed stationary gauge-field configuration, we derived an 
expression for the diquark contribution to the free energy that cannot be evaluated 
analitically but has a definite positive sign. 
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Integrating on the space of stationary gauge-field configurations the effective 
action in the saddle point approximation, by means of an expansion in powers of 
the gauge coupling constant, a gap equation is obtained compatible with previuos 
results. The gap is dominated by static chromomagnetic fields. 

We think, however, that a perturbative expansion in the gauge coupling con- 
stant cannot be fully justified even for large values of the chemical potential. 
Schematically, at the baryon density at which condensation of molecular diquarks 
is expected, the gauge couphng constant is presumably too large. For baryon 
densities for which the expansion might be justified, on the other hand, we ex- 
pect a BCS ground state, in which dibaryons have a size much bigger than the 
interquark separations, so that the dependence of the diquark structure functions 
on the spatial gauge fields cannot be ignored. 

We hope that our formulation should give a reasonable approximation to the 
QCD partition function for values of chemical potential of the order of the nucleon 
mass. Increasing the chemical potential we should meet chiral symmetry restora- 
tion, which according to our conjecture should be accompanied by the vanishing 
of the background field. If the phase transition is of first order, to determine its 
location we should compare the free energy evaluated in the present paper with 
that of the chirally symmetric phase in which the background field should vanish. 
But we notice that the latter cannot be simply obtained by setting the background 
field to zero in our equations. In fact in our saddle point approximation we disre- 
garded quasiparticles, appearing after the second Bogoliubov transformation, on 
the usual, reasonable assumption that they are separated from the vacuum by a 
large gap. If instead the background field vanishes, the particle antiparticle mix- 
ing in the action is again active, and if we construct the diquark field in terms of 
particles, we have no reason to expect that a gap exists for antiparticles. We must 
therefore proceed in a different way that we hope to illustrate in a future work. 
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A Superconductivity in many-body systems 

The phenomenon of superconductivity is explained in terms of Bose-Einstein 
condensation of fermion pairs, called Cooper pairs. The first suggestion in this 
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direction was made by Ogg |16], who observed persistent ring currents in solu- 
tions of alkali metals in liquid ammonia. But the importance his work was not 
understood and it did not have any influence in the development of the theory of 
superconductivity. 

The approach based on this idea is due to the Sidney group [2Z] , and it is formu- 
lated in the framework of the so called quasi-chemical equilibrium theory, in which 
there an equilibrium between formation and dissociation of Cooper pairs. The for- 
malism respects fermion-number conservation, but deals only with the ground state 
energy. In spite of its conceptual simplicity, calculations in this formalism are quite 
comphcated, and were completed only after the BCS theory [17] was published. 
We will also report this theory in the version of Bogoliubov and Valatin [MIES] , 
which does not respect fermion conservation, but introduces quasiparticles and is 
close to our approach from a technical point of view. 

A.l Quasi-chemical equilibrium theory 

The quasichemical-equilibrium theory is based on a variational calculation with 
a wave function constructed in terms of the wave function of the pair 

of the i-th and j-th among fermions 



*^ = ^E sgn(a)v.[7r(l),7r(2)]...^[7r(iV-l),7r(iV)] = C^(pf^) (A.l) 



Cat and a normalization constant, iSjv is the symmetric group of order and 
sgn((T) is the sign of the permutation a. The fermion Hamiltonian is 



with obvious meaning of the symbols. Since the Cooper pairs undergo Bose- 
Einstein condensation their total momentum is zero. Assuming that also the spin 
is zero, we get simple pairing: the fermions of the pair have opposite momenta 
and opposite spins, so that the Fourier transform (p of the pair wave function can 
be written function of a unique variable 



Here, and in the following, according to the standard notation, the sum over k 
implies also the sum over spins according to the previous identification. So that, 
if the probability of finding a fermion with momentum k is 




(A.2) 




(A.3) 



Pk 




,2 ' 



(A.4) 



24 



the condition of having fermions is 



J]p. = iV. (A.5) 

k 

To minimize the expectation value of the Hamiltonian under this condition, a 
chemical potential /x is introduced, so that one has to minimize the quantity 

£={Hn- fxNop) (A.6) 

where Nop is the fermion number operator. The above expectation value, whose 
evaluation EH is far from trivial is 



£ = p, + ^J2i(k,l\V\k,l) - {k,l\V\l,k)) pkPi 

k k,l 

+1 -k\V\l, -I) + O (j-) (A.7) 

k,i ^ ^ 



where 



e, := — , := • (A.8) 

Terms of order O (-^) have been neglected. They all respect fermion number con- 
servation, and can be disregarded below the critical temperature, but are essential 
in the determination of the critical properties of the transition to normal state. 

The first term is the contribution of the kinetic energy, the second of the 
density- density interaction, and the third of the interaction of fermions in one and 
the same pair. Because the density-density interaction is almost the same in the 
normal and in the superconducting state, to simplify the calculation it is accounted 
for by a renormalization of the chemical potential 

/^eff = - ^ 5^ ((A;, /| V|A;, - {k, l\V\l, k)) pi . (A.9) 
I 

Variation with respect to (p at constant fi^ff gives 

{ek - fics) Vk + ^{vl - I) = (A.IO) 

where 
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is the gap function. The solutions 



inserted in the definition of the gap function give the gap equation 



{k,-k\V\l,-l) 



I 2 



2 



A;. (A.13) 



In order to get a close solution some approximations are needed. First, a separable 
form is assumed for the matrix elements of the potential, that is 

otherwise. 

is the volume of the system. This is a crude but reasonable approximation 
of the electron-electron interaction due to phonon exchange. Second, since the 
contribution to the integral comes essentially from ~ jj^es ~ 2a7' "^^ere kp is 
the Fermi momentum, one can set ^ A^^ = A. Then the solution of the gap 
equation is 

^-2a;exp(^-^J • (A-IS) 

It is important to observe, in connection with our problem of including Cooper 
pairs and unpaired fermions at the same time, that, for \ek — fies\ > we assume 
Afc = 0, but need not to vanish. Actually in general ipk must not vanish for 



< /Ucff — in order to fulfill the condition f lA.Sp on fermion number. 



2m 

We conclude this subsection by discussing the effect of the coupling to a mag- 
netic field. The Cooper pair structure function depends on the applied gauge field, 
and the interaction can be computed in a perturbation series, that is 

= V7o + e <^i(A) + 0(e2) . (A.16) 

The first correction (fi (A) (e being the electric charge and A the vector potential) 
describes the Cooper pair magnetic susceptibility, and is essential in the expla- 
nation of the Meissner effect. Now the electromagnetic coupling is small with 
respect to the phonon coupling, which binds the electrons in a Cooper pair [37] , 
and therefore the perturbation expansion is justified. In the case of QCD it is, 
instead, the gauge interaction which binds quarks into diquarks, and therefore, 
in general, the dependence of the diquark structure function on the gauge fields 
cannot be neglected. 
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A. 2 Bogoliubov transformations 

The Bogoliubov transformations corresponding to simple pairing in standard 
notations are 



a 



- ^kC'k- VkC_k, ak = UkCk- VkC'_k (A. 17) 

where c|., Ck are creation- annihilation operators of the fermions in the system and 
the parameters m, u, not to be confused with the upper and lower spinor compo- 
nents, must satisfy the normalization conditions 



2 I 2 



1. 



Hi 



The transformed Hamiltonian is 

H' - = S + Hn + H20 
where Nop is the fermion number operator 

£ = J^(efc - fi)vl + ^ ^kUiVkVi {k, -k\V\l, -I) 



int 



kl 



k 



(cfc - n){ul - w^) - 2 ^ UkUiVkVi {k, -k\V\l, -I) 
(cfc - fi)uk Vk + l:{ul- vl) ^ uivi {k, -k\V\l, -I) 



O-k^^k 



(A.18) 
(A.19) 

(A.20) 
(A.21) 
(A.22) 



and the density- density interaction has been neglected. These terms have a close 
correspondence with those of the transformed action (14.151) . We have not written 
Hint- We only mention that it contains monomials of operators of power higher 
than 2 which do not conserve fermion number but are of order ( essential in the 
study of the phase transition to normal state). If we set 

fk 1 



Vk 



Uk 



(A.23) 



we see that the vacuum energy £ is identical to that found in the quasichemical 
equilibrium theory. At its minimum, H20 = 0, in perfect analogy with the results 
we got in the relativistic case, and 



(A.24) 
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so that the Bogohubov-Valatin method gives directly also the spectrum of quasi- 
particles which in the quasichemical equilibrium theory has to found separately. 
Introducing the parametrization (1A.23|) into the definitions (]A.17I) we can rec- 
ognize the form of our relativistic transformations. If we make the Bogoliubov 
transformation time-dependent, we can conserve fermion number by the help of 
compensating fields. The development of this approach for many-body systems 
can be found in ||48j. 



B Derivation of the diquark effective action 
B.l Pfaffians 

We first need to recall some basic facts about pfaffians. The interested reader 
can find a discussion about the properties of pfaffians and their relation to the 
Gaussian Berezin integrals in the detailed appendices of [19] , together with similar 
properties of determinants, permanents and hafnians. 

Let A = {Aij)f^^i be a be a 2m x 2m antisymmetric matrix. We define the 
pfaffian of A by 

= 5Z Sgn(a) A,(i),(2) • • ■ A,(2„^^1V(2H (B.l) 

where is the symmetric group of 2m elements, and sgn(cr) is the sign of the 
permutations a. Then 

(pfA)2 = detA (B.2) 

and for any 2m x 2m matrix X 

pf [XAX^) = (det X) (pf A) . (B.3) 

If A is invertible 

pf(A-^) = (pfA)-^ (B.4) 
Consider a partitioned matrix of the form 



M 



A 


B 




D 



(B.5) 



where A, B, D are matrices of sizes 2mx 2m, 2mx 2n and 2nx2n, respectively, with 
elements in a commutative ring with identity, and A and D are antisymmetric. If 
A is invertible, then 

pf M = {pf A) pf{D + B^A-^B). (B.6) 
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If D is invertible, then 

piM = {pfD)pf{A + BD-^B^). (B.7) 

Now let X15 • • • )Xn be the generators of a Grassmann algebra and A be an 
antisymmetric n x n matrix. Then the Gaussian Berezin integral provides a rep- 
resentation for the pfaffian: 

/" , , /I V"^ A \ I pf A if n is even ,^ ^, 

(2 g ) = [o it „ id odd . 

B.2 Evaluation of the transfer matrix 

In this section we refer to the fermion transfer matrix at non-zero chemical 
potential /i, in an arbitrary gauge, as defined in (13. 3p . with the notations given 
in (1321) and fl33|l . 

Let \uv) be the coherent state associated to the fermion operators u and i), 
that is 

\uv) = exp{—uu^ — t>£>^)|0) (B.9) 

where u and v are Grassmann variables. We shall make use of the completeness 
relation 

1= / dudu*dvdv* l^^^^^^l (B.IO) 



{uv\uv) 

by the help of the Berezin integration on Grassmann variables. 
The scalar product of two states is 

{uiVi\u2V2) = exp {ulu2 + vlv2) . (B-H) 

The evaluation of the matrix element of the transfer matrix between coherent 
states was already performed in [27] with the result 

{utVt\Tt,t+i\ut+iVt+i) = 

exp (u;nJv; + Vt+iNt+^Ut+i + u*f/o,te2^Mt+i + vlUl^e-^^Vt+i) . (B.12) 

Here we are interested in the evaluation of the matrix element 

X := (A,-^t|7;,m|A+i,-^m) • (B.13) 

Our procedure goes through the introduction of two complete sets of coherent 
states 

dutdu^dvtdv*dut+idu^_^_idvt+idv*^-^ 

—, , —{utVt\Tt,t+i\ut+iVt+i)—, . ^. (B.14) 

{UtVt\UtVt) {Ut+lVt+i\Ut+lVt+i) 
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Also to compute the matrix element {uiVi\T>^J^) we insert a complete set of co- 
herent states as follows 

{u,v,\V,:F)= [ du,du;dv,dv* (^i^i|exp(^V2)|^...)(^...|exp(^t)|o) 

J {U2V2\U2V2) 

According to our definitions (14. 8 p and (14. 2 p 

i)^ = u^R^V^R*h^ + vTR^V^R*^J^'^v - 2u^R^V^R*^T'^v (B.16) 
This means that 

{UiVi\exp{'Dy2)\u2V2) = {UiVi\u2V2) 

X exp (^ulR^^V^R*^2ul + ^V2J^R^V^R*^J^^V2 - u\R^V^ R*'^ J^'^ . (B.17) 
Similarly 

{U2V2\ exp(J^^)lO) = exp [ulF^l) . (B.18) 

Therefore 

{uiVilV, = exp ^-ulR^V^ R*^u{^ J du2du*2dv2dv2 exp {vi^F^v'^ 

xexp (\}2FR^V^R*^F'^V2 - u\R^V^R*^F'^V2 + u\u2 + vlv2 - u\u2 - vlv^ . 

(B.19) 

The integrations on the variables U2 and produce, respectively, the constraints 
u\ = ul and V2 = ulJ^^ , and we arrive at the result 

{uiVi\V,J^) = exp (^^ulR-^^V^{R-'^)^ul + ulJ^^vl^ . (B.20) 

By using all these intermediate steps we get 

X = y dutduldvtdvldut+idul^^dvt+idvl^^ 

fl _ 1 _i 
X exp ( -Ut{Rt'^)2'DtRt ^ut + VtJ'tUt - u*Ut - v*Vt 



X exp (u;n}v; + Vt+iNt+iUt+1 + u;Uo,te^''ut+i + v;U*te '^^Vt+i ^ 



(B.2i: 



30 



The integrations on the variables f j*^]^ and Vt produce, respectively, the constraints 
Vt+i = — and V* = —J^tUt, and, using the definition of J-jv,* given in f l4.26p . 
we arrive at the expression 



u 



t+l 



(B.22) 



As the next step we perform the Berezin integrations on Ut+i and u*^^ to get 



N,t+l) 



dutdu^exp ( -Ut{R^ ^)^VtRi ^Ut 



X exp 



X exp ( ^u;Uo,te^''iTl^t_,,] 



'r;^,vi,{r;^,)Hj^^ 



N,t+1/ 



. (B.23) 



These final integrals produce the pfaffian of a partitioned matrix and can be com- 
puted by using the expressions (IB. 61) or (]B.7|) to get 



X=det(J-^,+Jpf ^ 

X pf {Uo,te'^iTh,^,)-'R;M+iiRt+i)Hj'kt+^^^^ 

X [J^ N,t + UQ^ti^lf^t+iT^^t+iUl^t^ 



(B.24) 



where the last pfaffian can be re- written, by using formula flB.31) . as 



X = det(J'j 



■t 

N,t+1, 



pf 



{R;^)^VtRP det[J-^, + UoAJ'kt+ir'J'liUlJ't. 



xpi\[{R;^)mRt 



"21-1 



X 



[J'N. + UoAJ'kt+iy'J'liUlJ'tr^} ■ (B.25) 



Remark that det {J-'lf ^^^) = det(J^j^^^^f/Q J so that the product of this determinant 
with the other one appearing in flB.25p can be written as the determinant of the 
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product, which is exactly i?j+i,t according to f l4.27p . so that 



X=det(E,+i,i) pf [RfY^VtR] 



X pf 



+ e'^E;^,,R;;\vl,{R-^S^E;^,, 



(B.26) 



By using the relation fl4.34p and the formula fIB.Sp . we obtain, at the end, the 
expression 

X = det(E*+i,) pf (A) pf (D,-^ + e'^Q^i^{DliQt+i,t) ■ (B.27) 
We shall also need the normalization factor 

{V, J^lV, J^) = J dudu*dvdv* exp [^—u*u — v*v + u*T'^v* + vJ^u) 

which after the integration on v* and v and a rescaling of the variables becomes 
{V, T\V, T) = (det i?)"^ j dudu* exp (^u*u + ]^u*V^u* + ^mD^m^ (B.29) 



In conclusion 

exp(-S'feo) 



(det R)-^ pf (P) pf (P^^ + V^) . 



det(E,+i,) pf (P,) pf (I),-' + e^^Q^Vi AViQm,) 

t=o 

Lo/2-l 



(deti?),-ipf(I),)pf(Pr'+^^J 



JJ det(gt+i, 



t=o 



det(l + e^^Ag,;,Vi^l+iQm,* 



det(l + 



(B.30) 

(B.31) 
(B.32) 
(B.33) 



From this expression, as Sme is the part of S^o at P = and Sdq the rest, we easily 
derive that 



Lo/2-l 

Sme = " tr+ln(5t+i,t 

^ Lo/2-1 



t=0 



Sd, = 7: E ln(l + APj) - ln(l + e^'^AQ,7,ViAViQm,*) 



(B.34) 



(B.35) 
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